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Abstract. We present the long term evolution of the timing properties of the black hole candidate Cygnus X-l in the 0.002- 
128 Hz frequency range as monitored from 1998 to 2001 with the Rossi X-ray Timing Explorer (RXTE). For most of this 
period the source was in its hard state. The power spectral density (PSD) is well modeled as the sum of four Lorentzians, 
which describe distinct broad noise components. Before 1998 July, Cyg X-l was in a "quiet" hard state characterized primarily 
by the first three of these broad Lorentzians being dominant. Around 1998 May, this behavior changed: the total fractional 
rms amplitude decreased, the peak frequencies of the Lorentzians increased, the average time lag slightly increased, and the 
X-ray spectrum softened. The change in the timing parameters is mainly due to a strong decrease in the amplitude of the 
third Lorentzian. Since this event, an unusually large number of X-ray flares have been observed, which we classify as "failed 
state transitions". During these failed state transitions, the X-ray power spectrum changes to that of the intermediate state. 
Modeling this PSD with the four Lorentzians, we find that the first Lorentzian component is suppressed relative to the second 
and third Lorentzian during the state transitions. We also confirm our previous conclusion that the frequency-dependent time 
lags increase significantly in the 3.2-10 Hz band during these transitions. We confirm the interpretation of the flares as failed 
state transitions with observations from the 2001 January and 2001 October soft states. Both the behavior of the PSD and the 
X-ray lag suggest that some or all of the Lorentzian components are associated with the accretion disk corona responsible for 
the hard state spectrum. We discuss the physical interpretation of our results. 

Key words, stars: individual (Cyg X-l) - binaries: close - X-rays: stars 

1. Introduction nent which is in part reprocessed by the cooler disk. The 

hot plasma is also believed to be the reservoir for an outflow 
The existence of distinct states for the X-ray emission of black which emits synchrotron radiat ion at radio w avelengths (and 
hole binaries (BHCs) has been extensively documented for possibly also m the x . raySj see | M arkoff et aljEool and refer- 



many systems (see e.g.j Zhang et al.|1997aHWilms et al.|200U ences ^rzm). While the general observa tional proper ties of 
i Homan et alj | 200l|) . A geometrically thin, optically thick ac- B HCs a re no w well esta blished (see, e.g. Ivan der Klislfl99l: 



cretton disk is thought to be responsible for the soft, less van- l^^lEoorl iRemillardlEool . their unifi catio n in terms of 
able X-ray flux dominating the soft state. In the hard state, physical model is uncertain JZdziarskilll 9991: iNowaket alJ 
an additional hot plasma, the "accretion disk corona", is pos- Especially, there exists no consistent model for explain- 
tulated to produce a hard, highly variable emission compo- ing me complex short term varia bility on time scales faster than 
1000 s, although several phenomeno logical models have been 
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One of the possibilities to further improve current accre- 
tion models is to study the long term evolution of the canonical 
hard state, since in this state all the important physical com- 
ponents - disk, hot plasma, synchrotron outflow - as well as 
prominent short term variability can be observed. The stabil- 
ity of its characteristic parameters, however, has yet to be sys- 
tematically studied with high spectral and temporal resolution. 
We therefore have performed regular monitoring observations 
of Cygnus X-l with RXTE from 1998 to 2001. This persis- 
tent high mass X-ray binary (HMXB) is the best candidate for 
such a hard state campaign: it is bright (~0.6Crab) and spends 
most of its time in the hard state. In addition to the canonical 
hard state properties, the 2-12 keV X-ray flux and the 15 GHz 
radio flux of Cyg X-l display the 5.6 d orbital period and a 
150 d p eriod which might be due to precession of a disk/jet 
system llPoolev et alJll9 99: Brocksop p et alJll99 9a). The hard 
state shows weak and persistent radio emission that has been 
resolved as a type of steady jet (Stirling et al. 2001). The radio 
emission turns off during the soft state, as has been well docu- 
mented during the last extended s oft state which occurred from 
1996 May to 1996 August JZhang et alll997aUCui et alll997t 
ICui et all 1 9981: IWen et all200ll) . 

The hard state PSD of Cyg X-l has often been observed 
and modeled, the source serving as the BHC prototype. In 
addition to the flat-topped hard state noise, the occasional 
presence of low frequency noise on timescales longer than 
1000 s and of a probably associated weak quasi-periodi c os- 
cillation (QPO) at ~0.04 Hz has been reported jAngelini et alJ 
1 19941 IVikhlinin etall Il994l) . Early descriptions of the flat- 
topped PSD attributed the variability to Poisson distributed X- 
ray flares (shot noise) and generally approached the PS Ds as 
bein g due to a singl e broad-band proce ss jNolan et alJ ll98ll 
Lochneretal. 1991; Pottschmidt et al. 1998, see, however, 
Belloni & Hasinger 1990b). In recent years, however, it has be- 
come increasingly clear that the hard state PSDs of Cyg X-l 
and other BHCs sho w several "breaks", i. e ., frequencies wh ere 
the PS D steepens dNowak et al.1 Il999at iBelloni & Hasi nger 
1990a, and references therein), and that this complex structure 
has to be described by multiple components ( Smith et al. 1997; 
[ Berger & van der Klislll998t iNowaklEoOot iNowak et al.lbo02t 
IBelloni et alJl2002l) . Often, the ba sic component of such mod- 
els is still a shot noise component JPoutanenl2002tlSmith et ail 
1997). 

While on first sight the hard state PSD of Cyg X- 1 is domi - 
nated by two distinct broad components (Gilfan ov et all 19 99). 
strong indications were recently found that it can be described 
remarkably well with four broad noise components that yield 
the main contribution to the root mean square ( rms) amplitude 
in the 0.001 Hz to 100 Hz range JNowakKO OO). Occasionally, 
narrow features representing QPOs are required. A previously 
observed "QPO" peak, detected near 1 Hz with Q ~ 1 and 
~ 10% rms amplitude (Rutledge et alJll999l) . can possibly be 
identified with one of the above mentioned broad components 
being very prominent. During transitions to and from the soft 
state, broader QPOs (Q ~ 0.5-1, 3-17 % rms amplitude ) be- 
tween 3 and 13 Hz have been observed dCui et alJll99^) . As 
we find below, their parameters suggest that they correspond 
to the noise components in the hard state, with varying cen- 



tral frequencies and amplitudes. The normalization and flat- 
top break frequency - associated with the broad component 
that has the lowest characteri stic frequency - have long been 
known to be anti-correlated jBelloni & HasingeJll990bT) . In 
addition, the two lowest characteristic frequencies associated 
with broad compone nts have been observed to be correlated 
JGilfanov et alJll 999l). This latter correlation was discovered 
by Wiinands & van der Klis ( 1999) by comparing the flat top 
break frequency with the broad hard state "QPOs". 

A four-peak deconvolution of hard state black hole PSDs is 
especially interesting in the context of the question of whether 
any of these BHC noise features have the same origin as the 
generally much narrower and faster QPOs seen in neutron star 
binaries (e. g., the so-called horizonta l branch QPO and twin 
kHz OPOs:IWiinands & van der Klisll 19991: IPsaltis et alJll999t 
iNowaklfeoOOt Ivan Straaten et alJl2002l) . The aim of our PSD 
fits is therefore to identify the main noise components and to 
characterize their evolution in order to provide new data for 
constraining accretion models. 

In this work we concentrate on describing the variability 
properties of Cyg X-l in terms of the power spectral density 
(PSD), the Fourier-frequency dependent X-ray time lags, and 
simplified spectral models. Our aim is to characterize the accre- 
tion flow in the hard state. The description of the X-ray spectra 
in terms of more sophisticated accreti on disk corona model s 
will be published in a separate paper jGleissner et alJ l2003). 
This paper is organized as follows: In Sect. |3 we describe our 
observing campaign, the data screening, the computation of the 
power spectra, and the procedure of fitting multiple Lorentzian 
profiles to the power spectra. We then present our results in 
Sect. |5J starting with the long term evolution of the power 
spectrum (Sect. 13 - II . followed by a discussion of the failed 
state transitions (Sect. 13.21 1. the transition into the 2002 soft 
state (Sect. 13. 3> . and the behavior of the narrow Lorentzians 
(Sect. l3~4l . In Sect.@]we summarize our results. 

2. Observations and Data Analysis 

2.1. The RXTE Monitoring Observations of Cyg X-1 

The RXTE data presented here were obtained with the 
Proportional Counter Array (PCA; Jahoda et al. 1996) and 
with the All Sky Monitor (ASM; iLevine et all 1 19961) . us- 
ing the pointed RXTE observation programs P30157 (1998), 
P40099 (1999), and P50110 (2000-2001). We used the stan- 
dard RXTE data analysis software, FT00LS, Version 5.0, to cre- 
ate lightcurves with a time resolution of 2~ 6 s (P30157) or 2~ 8 s 
(P40099, P501 10) in multiple energy bands. The energy bands 
extracted for P30157 and P40099/P501 10 present the closest 
match possible for the available data modes (see Table [Q, and 
inclu de the bands used in our analysis of the X-ray time lags 
(Pottsch midt et alJ l2000. 2001). With a few exceptions noted 
below, the PSD analyses presented here were performed using 
summed lightcurves in the ~2keV to 13.1 keV energy range. 

The P30157 data set (weekly monitoring in 1998 with a 
nominal exposure time of 3 ks per observation) is the same 
that was used previous ly for the study of th e evolution of X- 
ray time lags (Pottschmidt et al. 2000, 2001) and the evolution 
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Table 1. Energy bands and data modes used. 



PSD 



lag/coherence 
low high 



PC A Epoch 3: data taken until 1999 March 22 

data mode B.16ms_46M_0.49_H 

channels 0-35 0-10 23-35 

energy [keV] -2-13.0 -2-4.2 8.3-13.0 

PC A Epoch 4: data taken after 1999 March 22 

data mode B.2ms.8B_0_35.Q 

channels 0-30 0-10 20-30 

energy [keV] -2-13.1 -2-4.6 8.4-13.1 



of X-ray spectra ( Smit h et al]l2002l) . The data set P40099 (bi- 
weekly monitoring in 1999 with a nominal exposure time of 
lOks per observation) presented here has also b een included 
in the time lag study of Pottschmidt et al. (2001). The P501 10 
data (biweekly monitoring with a nominal exposure time of 
20 ks per observation) have not been published before. Starting 
in 1999, simultaneous observations at radio wavelengths per- 
formed with the Ryle telescope in Cambridge, U.K., are also 
part of the campaign. In addition, daily radio fluxes were ob- 
tained with the Ryle telescope. 

For all observations, high resolution lightcurves and X- 
ray spectra were generated using standard screening criteria: a 
pointing offset of < 0°01, a SAA exclusion time of 30 minutes 1 , 
and a source-elevation of > 10°. At times of low source lumi- 
nosity, when the "electron ratio", i.e., a measure for the particle 
background in the PCA, was not influenced by the source itself, 
we also imposed a maximum "electron ratio" of 0.1. Good- 
time intervals (GTIs) were then created for all available PCU 
combinations, and from those suitable GTIs were selected for 
extracting high time-resolution data for the data modes and en- 
ergy bands given in Table ^ F° r each GTI file, we generated 
the background spectrum using the FTOOL pcabackest and 
the "sky-VLE" background model. The background in a given 
energy band was then taken into account in determining the 
PSD normalization. 

As a result of the data screening, ~2ks of usable data 
were left for each of the P30157 observations and ~7ks for 
P40099 and P50110. Note that due to RXTE scheduling con- 
straints the full nominal exposure time was only rarely reached. 
Observation dates and exposure times are included in Table|2] 
The table caption given here contains a description of all pa- 
rameters, the full table is available in electronic form at EDP 
Sciences. 

2.2. Computation of the Power Spectra 

The computation of the po wer spectra for all energy bands 
given in Table □ follows iNowak et alJ Jl999ah and iNowakl 
<200(J. 

For the P30157 data set, PSDs in the 0.002 to 32 Hz range 
were computed, while the P40099 and P501 10 PSDs reach up 
to 128 Hz. This discrepancy arises from the different maximum 



time resolutions of the binned data modes that are available 
for these data sets (see Table[Q. We decided to use the binned 
modes with their comparatively moderate time resolutions, be- 
cause they allow us to study several different energy bands 
(see Sect.[^} which cannot be done with the "single bit" data 
modes (250 //s time resolution) that are also available. 

All PSDs used in this paper are pr ese nted in the nor- 
maliza tion of lBelloni & Hasingerl i ll 99 0b) and Miya moto et alJ 
Jl992l) . where the PSD integrated over a certain frequency in- 
terval equals the square of the relative contribution of that fre- 
quency interval to the total rms noise of the source. Where 
possible, we adjusted the normalization to the background 
corrected count rate. The observational deadtime corrected 
Pois son noise was subtract ed before t he no rmalization. See 
Vikhlinin et al J I I 19941 and Zhang et alJ ll 19951) f or a discussion 
of the dea dtime influence on th e PSD in general. ljernigan et alJ 
J2000h andlZhang et al] ill 9961) for the case of the PCA detector, 
and iRevnivtsev et al.1 1 120001) for an application to PCA mea- 
surements of Cyg X-l above 100 Hz. For mos t observations 
we use d the first two components of Eq. (6) of iTernigan et al.l 
(2000), namely the approximati on of the "general paralyz- 
able deadtime influence" (Zhang et al 1199.4 Eq. (24)) plus the 
deadtime caused by "Very Large Events" (characterized by the 
PCA VLE count rate). For frequencies < 100 Hz the correc- 
tion of the noise level due to these two components amount 
to only a few percent. Higher order corr ections or a fit of these 
compone nts to the PSDs a s applied by Jernig an et alJ (J2000) 
or IRevnivtsev et alJ J200 0) are therefore not necessary here. 
We note that the discussed features in the PSDs have rms am- 
plitudes well above the typical deadtime str uctures (see, e.g., 
1 l2000l or Fig. 4 of IRevnivtsev et all 



Fig. 5 of Jernigan et al. 



For a source as bright as Cyg X-l this is a very conservative crite- 



A slightly different PSD normalization and deadtime cor- 
rection has been applied to all observations performed after 
2000 May 12, when the epoch of enhanced PCU background 
count rates started. Here, a good background model was not yet 
available and thus we could not background correct the PSDs. 
This correction is only on the order of a few percent and does 
not influence appreciably our results (we note that this problem 
is, h owever, a serious p roblem for sources with a low count 
rate; iKalemci et ai1l200ll) . Furthermore, as we do not have ac- 
cess to the VLE count rates from individual PCUs, it becomes 
impossible to correct the PSD for the VLE specific deadtime. 
Therefore, we only applied the general paralyzable deadtime 
correction. 

For each observation we created three PSDs by averag- 
ing the individual PSDs from lightcurve segments with dura- 
tions of 512 s, 128 s, and 32 s. Logarithmic frequency rebin- 
ning was performed on the same frequency grid for each of 
the three resulting PSDs. This grid ranges from 2~ 9 to 2 5 (or 
2 7 )Hz (-0.002 to 32 [or 128] Hz), with different df/f values 
for different frequency ranges: 25 logarithmically spaced bins 
were used from 0.001 Hz-0.1 Hz, 50 bins from 0.1-20 Hz, and 
5(15) bins from 20-32(128) Hz. The final PSD was then cre- 
ated by selecting for each grid bin the highest signal to noise 
value among the three PSDs. 
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Fig.l. a) Best fit (solid line) to a typical 0.002-128 Hz hard 
state power spectrum of Cyg X-l (P40099/18, 14off, 1999 
September 13, diamonds) using multiple Lorentzians. The fit 
was performed for the energy range from ~2 to 13.1 keV. It in- 
cludes four broad (dashed lines) plus two narrow (dash-dotted 
lines) Lorentzian profiles and a power law component (dotted 
line), b) The same as above but here the PSD has been mul- 
tiplied by frequency to illustrate the peak frequencies of the 
broad noise components (see Sect.@]for a discussion of the im- 
portance of these frequencies), c) Residuals between data and 
model in units of the standard deviation of the PSD data. 



2.3. Modeling the Power Spectra 

As was shown bv lNowakl (2000). multiple Lorentzian profiles 
can provide a good description of the typical hard state vari- 
ability of Cyg X-l and GX 339-4 in the 0.001-100 Hz range. 
Inspired by this success, we applied multiple Lorentzian pro- 
files to model the PSDs from our monitoring campaign and 
generally obtained good fits (see Tableland the typical exam- 
ples displayed in Fig.^and Fig.|5Jl. Note that power spectra ob- 
tained during the soft state, e.g., during the 1996 soft state, are 
not well fit with the multi-Lorentzian model (see also Sect. 13.21 
below). Our fitting approach uses standard x 1 minimization 
with the uncertainty of th e free fit parameters de termined using 
the usual prescription of Lampt on et al.l 41976). Unless noted 
otherwise, we give the uncertainty in terms of lcr errors. 




Frequency [Hz] 



Fig. 2. a-c) Another typical fit, in this case to one of the P30157 
observations (No. 11, 1998 February 20) in the 0.002-32Hz 
frequency range. The best fit consists of four broad and one 
narrow Lorentzian profiles. 



We describe the power spectra as the sum of Lorentzian 
profiles of the form 



Uf) = n- 



(1) 



where /j is the resonance frequency of the Lorentzian, ~ 
fil A//,fwhm its quality factor 2 , and Rj its normalization con- 
stant. Note that our definition of Q is consistent w ith the com- 
mon use of the O- value (e.g.. | Belloni et alJl2002t). but differs 
from that used by iNowaketalJ Jl999al) and lNowakl (|2000) 
where Q = /i/A/ i(HW HM with A/J-hwhm = was 
used. 

Integrating this profile over the frequency range from zero 
to infinity gives its total rms amplitude: 



'1 tan-'(-2&) 
rms, = R t | 

i 2 71 



1/2 



(2) 



A further important quantity of the Lorentzian is its peak fre- 
quency, i.e., the frequency where its contribution to the total 



" A/|_fwhm is the full frequency width of the Lorentzian at half of 
its maximum value. 
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rms variability is at its maximum: 



Vi = fi 



1 + 



4G? 



1/2 



(3) 



As it was recently found that v, seems to be the important 
parameter in terms of frequency correlations INowakl 12 000: 
NowaketalJ Eool Ivan Straaten et all Eool iBelloni etaTI 



v'ak ( 



2002), we will use v,, and not /;, to describe the location of 
each Lorentzian, unless specifically noted otherwise. 

Alternative models for the lowest frequency Lorentzian, 
such as using a zero centered Lorentzian, i.e., a sh ot noise com- 
ponent featuring a single relaxation time scale ( Be lloni et alJ 
120021) . may in principle provide fits of comparable quality (note 
that in the case of the observation analyzed by iNowa kl l2000l 
a zero centered Lorentzian provided a significantly worse fit). 
The multi-Lorentzian fits, however, provide the most uniform 
description as well as the possibility to compare the noise com- 
ponents with those found in neutron star binaries (see Sect.|4}. 

In almost all of our PSDs, four broad Lorentzians are ap- 
parent (Figs.n an dE) which we call L\, Lq, L3, and L4. The 
fourth broad component, L4, is usually present a bove 32 Hz 
and ha s also been noted by Revnivtsev et al.l J2000l) and lNowakl 
(2000). This component is also seen in other black hole or neu- 
tron star sources (van Straaten et al. 2002; Belloni et al. 2002). 
In the P30157 data set only its low frequency part can be seen. 
However, with the knowledge from the P40099 and P50110 
data sets, we constrained its shape, especially its approximate 
width, and included this component in the P30157 models. We 
will discuss the issue of possible systematic errors resulting 
from this approach later. 

We generally achieved a good fit using four broad 
Lorentzian components, L1-L4, with typical peak frequencies 
of ~0.2 Hz, ~2 Hz, ~6 Hz, and ~40 Hz. In our first modeling 
attempts we left the Q-factors of the Lorentzians as free param- 
eters. They were found to vary considerably from observation 
to observation: Q\, Q2, and Qj ranged from ~0. 1-2.8, Q\ was 
even less constrained. However, later experiments showed that 
fixing each Q-valne to a representative value for each broad 
component does not significantly alter the fit results. This can 
be expected, since the broad components overlap over large fre- 
quency intervals (Figs.^5 and|5J)). We therefore set Q\ = 0.25, 
Q 2 = 0.3, Q 3 = 0.3, and Q 4 = 0.5 for all fits. This resulted in a 
good description of most PSDs. There are only 22 cases out of 
the total 130 PSDs where the fit significantly improved when 
one or more of the Q-factors were treated as a free parameter. 
In these cases, which mainly coincided with the "failed state 
transitions" further discussed in Sect. !3.2l below. we accepted 
the model with a variable Q as the best fit. For the low signal 
to noise energy resolved PSDs of the P30157 data we fixed Q4 
and /4 to their best-fit value of the non energy resolved PSD. 
In addition to L\ through L4, a power law of the form 



PL(/)=/? pI /- 



(4) 



is required in ~40% of the observations. The majority of 
these power laws are required to reduce the residuals below 
~0.01 Hz. They have a typical slope of ~1.3 and a normaliza- 
tion of ~10~ 4 . In some instances the power law is the main 



contributor to the total rms amplitude. As we show in Sect. 13. 21 
these observations mainly coincide with "failed state transi- 
tions". Generally, these observations are the same where we 
also had to leave the Q-values of the Lorentzians a free pa- 
rameter. This can be expected as the contribution of the power 
law at higher frequencies slightly contaminates the other fit pa- 
rameters. In other observations - normally those close to the 
"failed state transitions" - a clear deviation from the fit with 
broad Lorentzians is also present at low frequencies. Due to 
the logarithmic binning of our data, in these cases the devia- 
tions were quite often only seen in the lowest few frequency 
bins and thus there was not sufficient information to reliably 
constrain the shape of the deviation. For consistency with the 
fits performed for the "failed state transitions", we chose to de- 
scribe these observation with the same PSD model that we also 
used for the intermediate or soft state, thus implying that the 
reason for the low-frequency excess in these observations is the 
same as that for the excess during the "failed state transitions". 
Since the uncertainty of the power law slope is rather large or 
undeterminable due to numerical reasons for these latter obser- 
vations, we do not give specific uncertainties. We also caution 
that the difficulty in constraining the power law also results in 
an anti-correlation between a p \ and the power law normaliza- 
tion R p i, which we view as being largely systematic in nature. 

In agreement wi fh lNowakl ( l2000l) . lNowak et alJ ( l2002l) . and 
Bell oni et alJ II2002I) . we often see additional subtle but statisti- 
cally significant substructures. Sometimes these structures are 
sharp (Figs.n an d|3 or weak and broad. In this case, further 
Lorentzian components (L a dd,;) were added to the basic PSD 
continuum. We generally retained these components if their ad- 
dition resulted in a significant improvement of the best fit^ 2 , 
or if they were very clearly present in the PSDs 3 . Residuals 
at frequencies slightly higher than the peak frequency of L\ 
tended to be narrow. We modeled these features with one or 
two narrow Lorentzians, L a dd,i and L a dd,2 (see Fig.^and Fig. [2}. 
Typical peak frequencies for these QPOs range from 0.1 to 
0.6 Hz, and one might speculate about their relationship to the 
horiz ontal branch o scillation peaks in neutron star X-ray bina- 
ries ([Wiinands & van derKlislll999l Table 4). In addition, be- 
low L\ a broad, weak (< 2cr) residual at ~0.03 Hz was often 
seen, which we designate L a dd.3- The total contribution of L a dd,i, 
i a dd,2, and L a dd,3 to the rms is on the order of a few percent at 
most. See Sect. 13.41 for a discussion of the behavior of these 
additional components. 

Using this modeling approach, we obtained values for the 
reduced x 1 m the range of 1-2. The best-fit parameters of all 



3 As mentioned in Sect. 12. 21 the PSDs are logarithmically rebinned 
to give a good description of the continuum power. This rebinning 
often had the result that the narrow features were confined to only one 
or two frequency bins. As a result, their modeling with a Lorentzian 
could not improve the x 1 although the feature is clearly present in 
the unbinned power spectrum. To be at least able to note the presence 
of such a feature, we nevertheless included them in the fit; however, 
only their frequency and power can be constrained, while their width, 
as represented by the g-value, can only be described as "narrow". 
These features will generally have Q > 50 in Table |2| Because of 
these complications, we also do not give any formal uncertainty for 
the narrow Lorentzians in Table|5| 
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Table 2. Results of fitting multiple Lorentzians to selected power spectra of the Cyg X-l RXTE/PCA monitoring observations 
from 1998-2001 in the frequency range from (1/5 12) Hz to 32 Hz (P30157) or 128 Hz (others), for the 2-13 keV band (TableQ). 
The table is available in electronic form via http://www.edpsciences.org Apart from the fit results, it contains the following 
information for each observation: running number, proposal name (e.g., "P40099"), sub-ID within the proposal, PCUs turned off, 
date of the observation, and the exposure time in ks. The fit results for each observation are organized as follows: for each of the 
possible four broad and three thin Lorentzians the norm (/?), central frequency (/), and quality factor (Q) with their corresponding 
uncertainties are listed. In addition the power law normalization (A) and slope (a) are given if present in the model. Finally, the 
^ 2 /dof and;^ 2 ed values of the fit are listed. 



130 observations are given in electronic form in Table|2] Most 
of the remaining systematic residuals in these fits are probably 
due to the fact that the broad noise components have a more 
complex shape than a simple Lorentzian profile and/or that sub- 
harmonics are present. This is especially true for L2, which is 
often associated with <2cr residuals in the ~ 1-2 Hz range (see 
also Nowak 2000). The broad low-frequency residual L a dd,3 and 
other residuals at frequencies below v\ might also be due to this 
effect. A further source of weak systematic residuals might be 
our choice of fixed Q-factors. 

After obtaining the final best fit we computed the total rms 
variability amplitude relative to the mean count rate by sum- 
ming the contributions of all power spectral fit components. 
The best fit Lorentzian profiles were integrated from zero to in- 
finity (Eq. (|2l), while the power law contribution was obtained 
by integrating the power law over the frequency range from 
10~ 3 Hz to 200 Hz. Because of the different time resolution of 
the extracted lightcurves this approach is more systematic than 
the direct measurement of the rms variability from the X-ray 
data. 

2.4. Time Lags, Coherence, and Spectral Modeling 

Although the emphasis of this paper is on the evolution of the 
short term timing properties of Cyg X-l in terms of the PSD, it 
is obviously important to also characterize the source in terms 
of other quantities. 

For black hole candidates, the X-ray lightcurves in two 
energy bands are similar to each other. However, the X-ray 
lightcurve in a higher energy band generally lags that mea- 
sured in a lower energy band. This X-ray time lag depends 
on the Fourier frequency. A measure for the similarity of the 
X-ray lightcurves in these two bands is given by their coher- 
ence, which is again Fourier frequency depend ent. We com- 
pute both quantities using the formulae given bv lNowak et alJ 
Jl999al) . As we have shown previously, an especially interest- 
ing quantity is the average X-ray time lag in the 3.2-10 Hz band 
(Pottschmidt et al. 2000), which we will use to quantify the X- 
ray lag. 

To obtain a description of the X-ray photon spectrum, we 
describe each source spectrum as the sum of a power law spec- 
trum with photon index T and a m ulti-temperature disk-black 
body after Makish ima et alJ dl986l) . characterized by the tem- 
perature at the inner edge of the accr etion disk, kT\„. To this 
contin uum, a reflection spectrum after Magdziarz & Zdziarski 
il 199.4 is added. For this paper we concentrate on the X- 
ray spectrum below ~20keV by ignoring the HEXTE data. 
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Fig. 4. a) RXTE-PCA count rate (2.5-20 keV) for each point- 
ing, as well as b) RXTE-ASM count rate and c) 15 GHz radio 
flux, both binned to a resolution of 1 d to show the variability 
on timescales that are smaller than the orbital timescale. 

Such an empirical model is sufficient to roughly describe the 
most important spectral parameters. A full study of the evo- 
lution of the broad-band X-ray spectrum that includes the 
HEXTE data and models th e source spectrum in terms of th e 
Co mptonization models o f iDove et alJ Jl997l) . ICoppil {£999), 
and lPoutanen & Svensson ( 1996) will be presented in a forth- 
coming paper (Gleissne r et alJ2003l) . 

3. Evolution of the Power Spectrum 

In this section we will discuss the overall behavior of the PSD 
of Cyg X-l since 1998. The evolution of the best fit PSD pa- 
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Fig. 3. Spectral and timing parameters of Cyg X-l between 1997 and 2001. a) Photon index, b) X-ray time lag between 2^1 keV 
and 8-13 keV, averaged over the 3.2-10Hz frequency band, c) X-ray coherence between 2-4 keV and 8-13 keV, averaged over 
the 3.2-10Hz frequency band, d) Relative rms amplitude for the total 2-13.1 keV power spectrum, e) Relative rms amplitude 
and f) frequency of each of the four broad noise components. Crosses: rmsi, diamonds: rms2, circles: rms3, triangles: rms4. g) 
Daily averaged 15 GHz radio flux, determined with the Ryle telescope, h): RXTE ASM 2-12 keV count rate. The dotted vertical 
lines mark special observations discussed in Sect. 13.21 and Sect. 13.31 The radio flux and the ASM rate have been rebinned to a 
resolution of ~5.6 d, the orbital period of Cyg X-l, for clarity. 
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Table 3. Average values and standard deviations of the rms amplitudes, rms,, the peak frequencies, v,-, and the frequency ratios, 
Vj/vj, for the four broad noise components (i, j = 1, ... ,4). The values are derived for different time episodes surrounding the 
1998 May change and the 1998 July failed state transition. For the P30157 data the highest peak frequency, V4, is biased (see 
text) and is therefore not quoted in this table. 
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Fig. 5. Constancy of the peak frequency ratios over the mon- 
itoring campaign, a) Diamonds: V2/V1, stars: V3/V1, triangles: 
V4/V1. b) Stars: V3/V2, triangles: V4/V2, crosses: V4/V3. 



rameters, the X-ray time lag, and the X-ray coherence func- 
tion, as well as the photon index, T, from the spectral fits, the 
15 GHz radio flux, and the RXTE ASM count rate are shown 
in Fig. [5] Both the radio flux and the ASM count rate have been 
rebinned to a resolution corresponding to the orbital period of 
the HDE 226868/CvgX-l system (P nr b = 5.599829(16) d, see 
iBrocksopp et al.ll999bHLaSala et al.ll 998). This rebinning re- 



moves the effect o f the well known orbital variation of the X- 
ray and radio flux JPoolev et alll999HBrocksopp et all 19 99a: 
Wen et al. 2001 ) such that long term changes become more ap- 
parent. Due to the rebinning, however, many short flares that 
are not due to the orbital modulation are smoothed out. We 
therefore present the radio and ASM data again in Fig.@]with a 
higher time resolution, together with the RXTE PCA count rate 
determined from the monitoring observations. Not unexpect- 
edly, the ASM and PCA rates are well correlated, although their 
correlation is not perfect due to their different energy bands and 
the X-ray spectral variability. 

3.1. The Two Hard States 

The most striking result of our fits is the bimodal behavior of 
the total rms amplitude (Fig.|3Jl): Starting in 1998 mid- April, 
the overall variability of Cyg X-l on short timescales showed 
a decrease of ~7%, from an average rms amplitude of 36±1% 
to 29+1%. At the same time, the power law index of the X-ray 
spectrum slightly softened from ~ 1.75 to ~ 1.85, the charac- 
teristic frequencies of L\ through L3 shifted to higher values, 
and the X-ray time lag increased. Furthermore, a larger frac- 
tion of observations requires the power law component below 
~0.01 Hz. Apart from further "X-ray flares", which we will de- 
scribe in Sect. 13. 21 Cyg X-l has remained in this softer and less 
variable (on short time-scales) X-ray state. Table[3]summarizes 
the main properties of the PSD before, during, and after the 
1998 changes. 

The parameters for the X-ray spectrum and timing behavior 
measured before 1998 May are very similar to the canonical 
hard state values for this source. Our campaign shows that these 
values have not been reached since then, i.e., since 1998 May, 
Cyg X-l is exhibiting a behavior that is atypical compared to 
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Fig. 6. Normalization of the four broad noise components, /?,-, 
versus their peak frequencies, v, . The different plot symbols are 
attributed to the same components as in Fig. [21 

previously observed hard states of this source. In the following 
we describe these changes in further detail. 

The change of the total rms can be mainly attributed to the 
behavior of L3 (Fig. |3£): while the rms values of L\ and Lq 
stay roughly constant, L3 declines from its pre-May strength 
of ~ 17% to 10+3%, i.e., almost half of its original strength 4 . 
Especially during the X-ray flare in 1998 July, L3 cannot be 
detected at all. Since these events, L3 has mainly been present 
at a level between ~5% and ~10% rms, and was sometimes not 
detected at all. Typical examples for pre- and post- 1998 May 
power spectra are shown in Figs. [5] and ^ respectively. These 
figures clearly show the decrease in strength of L3, especially in 
the panels where the PSD has been multiplied with frequency. 

In addition to this large change in strength, the whole PSD 
shifts to higher frequencies (Table [3}. The two lowest peak 
frequencies, v\ and V2 increase by a factor of ~3 during the 
change, while V3 increases by a factor of ~2. This behavior is 
not mirrored by the peak frequency of L\. The sharp drop from 
about 60 Hz to 40 Hz at the end of P30157 in v 4 is due to the 
biasing of its value caused by the low frequency cutoff of our 
data. We note that although the peak frequencies change, their 
ratios are remarkably constant over the whole campaign - in 
fact, the 1998 May event is only barely discernible in the evo- 
lution of the frequency ratios (Fig. [5j». 

The nature of the change in the hard state timing behavior 
can be clearly seen when the variations in the variability am- 

4 Due to the high frequency cutoff of the P30157 data, the parame- 
ters of L4 are difficult to constrain during 1998. We found no indication 
for a bimodal behavior of this component, however, the large uncer- 
tainty of the fit parameters does not allow us to make any definitive 
statements. 



plitudes and peak frequencies are considered together. Fig. [6] 
shows the dependence of the normalization constant, Ri, of 
each broad component on the associated peak frequencies, v,, 
for our data. Regions in Fig. [6] where the data point density 
is higher are due to the hard state behavior before and after 
1998 April. 

Several interesting features are apparent in Fig. [6] The 
three lowest frequencies, vi through V3, form a continuum from 
about 5 10~ 2 Hz to 20 Hz. For all three Lorentzians, the nor- 
malization of the Lorentzian correlates with its peak frequency 
- higher v, also imply a lower normalization. Furthermore, 
in those frequency regions, where several Lorentzians can be 
found, e.g., around 0.5 Hz for v\ and V2, and around 5 Hz for 
V2 and V3, their normalizations seem to be the same. We stress 
that this is not due to a misidentification of the Lorentzians 
as these are always clearly distinguishable in the observations 
presented here. There are two distinct regions of the Z?,-Vj cor- 
relation: below ~ 3 Hz, where the figure is dominated by L\ 
and La, and above 3 Hz, where L3 dominates. Clearly, V3 and 
/?3 are much more strongly correlated than the parameters for 
the low frequency Lorentzians. We note that the sharpness of 
this correlation might in principle have been influenced by the 
32 Hz maximum frequency of our 1998 data set. However, the 
strong 7?3-V3 correlation is also present when only considering 
the P40099 and P50110 data. We speculat e that these correla- 
tions ar e similar to the effect first seen by iBelloni & Hasing erl 
(1990b) in terms of their broken power law analysis of the 
PSDs of Cyg X-l, where the break frequency of the PSD is cor- 
related with the total source rms variability. In our Lorentzian 
decomposition, v\ is roughly equivalent to the break frequency 
and shows a similar behavior. What is new, and not yet under- 
stood, however, is the behavior of L3. 

Finally, the highest peak frequency component, L4 (trian- 
gles), is clearly distinct from the lower frequency features and 
shows a different behavior. Its peak frequency and normaliza- 
tion are comparatively stable. The clustering of the data at 40 
and 60 Hz again reflects the fact that V4 has been systematically 
over-estimated during P30157. During P40099 and later, how- 
ever, no systematic changes are seen in V4, although our data 
should have been sensitive to such changes. 

We now discuss the properties of the Lorentzians in terms 
of other measured quantities (Fig.[7J. The distinct behavior of 
the Lorentzians when plotted against the X-ray photon spec- 
tral index, F, and the X-ray time lag gives us further confi- 
dence in the individual identifications of the Lorentzians. Both 
the X-ray time lag and the peak frequencies v\ through V3 
are clearly correlated with F: softer photon spectra imply a 
shift of the PSD towards higher frequencies and also higher 
X-ray time lags (Fig. [7J) and c). There is also a clear corre- 
lation between the average time lag and the Lorentzian peak 
frequencies (Fig. [7^). Although a correlation between photon 
index and fre quency has been noted before with more lim- 
ited data sets jdiMatteo & Psaltislll999t iGilfanov et alJll999t 
iRevnivtsev et alJ20 00: Now ak et all2002l) . this is the first time 
that a simultaneous correlation with the time lags has been 

found for Cyg X-l. A similar behavior has be en seen in a 

J 1 

more limited hard state data set of GX 339-4 (Nowak et al. 
2002). Recently, a comparative study of the "QPO-r" correla- 
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tions has been performed for several BHCs with narrow QPOs 
bvlVi gnarca et al.N2003l) . While the direct identification of one 
of our correlations of Fig.^b with the unified track displayed 
in their work is not easily possible (different spectral models 
were used to derive T), the parameter space is comparable. 
We do not see any turnoff, i.e., no decreasing photon index 
at high frequencies which seems to be a common f eature in 
some of the sources studied by Vign arca et alJ ( 120031) (mainly 
GRS 1915+105, GRO J1655-40, and XTE J1550-564). This 
is consistent with the picture suggested by those authors that 
the turnoff might be related to a transition between the hard and 
the very high state, since the latter is not observed in Cyg X-l 
(at least not in its canonical form, note, however, that the hard 
spectral component is generally present also in the soft state). 
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Fig. 8. Peak frequencies of L\, L2, and L3 as determined for the 
energy band from ~ 2^4.6 keV versus the peak frequency for 
the energy band from 8.4-13. Ike V. The symbols correspond 
to those used in Fig. [3] 



We note that there is no strong energy dependency of the 
power spectrum shape during the normal hard state phases. As 
an example, Fig.[8]shows the correlation between the peak fre- 
quency of L\, L2, and L3 for the two energy bands for the 1998 
data. To within their error bars, the peak frequencies of L\ and 
Li are identical for these bands. The peak frequency of L3 is 
higher by about 30% in the hard energy band, a trend that con- 
tinues when analyzing lightcurves in even harder energy bands. 
This trend is similar to that seen in t he Lore ntzian decomposi- 
tion of the PSD of XTE J 1650-500 JKalemci et all2003l) . 

3.2. Failed State Transitions 

In this section we study the behavior of Cyg X-l during the 
"flares" apparent in the ASM lightcurve of Fig. [3] During these 
flares, the rms amplitude decreases, the strength of the power 
law contribution relative to the Lorentzians increases, and the 
X-ray spectrum softens. Generally, the X-ray time lag also 
shows values that are much higher than those seen before and 
after the flare, and the coherence function drops. PSDs mea- 
sured during several of these events are shown in Fig.|9] These 
examples are significantly different from the "standard PSD 
shape" for the hard state as defined in Figs.^and[2l 

Crucial to the interpretation of these flares is the evolution 
of the PSD and the other timing quantities over the flares. We 
will concentrate on the flares best sampled by our monitoring 
observations. These flares, which were observed in 1998 July, 
1999 December, and 2000 November, are identified by dashed 
lines in Fig. [3] Additional smaller events were also seen, e.g., 
in 1999 September, in 2000 March, and in 2000 December. For 
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Fig. 9. a-d) Examples of power spectra that deviate significantly from the typical hard state shape as defined by Figs.[T]and|2 
The best fit broad noise components are distinguished by their line styles: dashed: L\, solid: L2 and L3, dash-dotted: L4. The 
remaining components and the total best fit are dotted. The observation dates are given in the plot, they are also marked by dotted 
lines in Figs.[3]and|4] 



these latter events, however, our observations did not sample 
the change of the PSD in sufficient detail. 

A typical example for the evolution during the first stage 
of these flares was seen in 1998 July (Figs. l9l and I10>: Here, 
the characteristic frequencies of the Lorentzians shift to higher 
frequencies and their relative strength changes. In contrast to 
the standard PSD, the L3 component is weak or missing and 
the peak frequencies of L\ and L2 are significantly enhanced. 
Furthermore, the X-ray time lag in the 3.2-10 Hz band in- 
creases (Fig.|3j), see also Pottsch midt et all2000l) . On the reso- 
lution of our monitoring, the transition back into the hard state 
mirrors that of the transition into the flare. 

We briefly note that during the 1998 July event there is a 
change in the energy dependency of the rms variability am- 
plitude (Fig. II It . Before the flare, the rms variability ampli- 
tude of L3 increased with photon energy, that is the PSD was 
flatter at higher energy bands as is typical for the hard state 
jNowak et alJl999al) . During the flare, when the spectrum was 
soft, the rms variability amplitude decreased with energy. This 
is best illustrated in the behavior of L\, where the total power 
in this component before the 1998 July event has a rather small 
energy dependence (Fig. m solid lines), while it is strongly 
energy dependent during the flare (Fig.^ dashed lines). 



During some of the flares the PSD evolution starts similar to 
1998 July but then continues until a strong low-freqency power 
law dominates the observed rms, instead of the Lorentzians. 
Typical examples for such a behavior are observations made 
on 1999 December 05 and on 2000 November 03, at the 
peak of their respective outbursts (Figs. 110b and c). Of the 
Lorentzians, Lo is strong and L3 is clearly present, albeit at 
the weak levels typical of PSDs observed after 1998 May. On 
the other hand, L\ is barely identifiable. Both observations also 
have quite large X-ray lags (e.g., ~8ms for 2000 November) 
and decreased coherence (e.g., 0.6 on 1999 December 05). 
The neighboring observations, however, display the charac- 
teristic double -peaked PSD similar to 1998 July, coherence 
close to unity, and lower lags of ~3 ms. The evolution of 
the Lorentzians through the 2000 November flare is shown 
in Fig. El Both, the PSD shape as well as the X-ray spec- 
tral parameters seen here ar e very similar to ea rlier observa- 
tions of Cyg X-l reported bv lBelloni et al.Nl996l) . Comparing 
the timing pr operties and spectrum of GS 1124-683 (N ova 
Muscae 19 91:lEbisawa et alJll994 [Miyamoto et al.lll994l) and 
GX 339-4 JMendez & van der Klisll997t) to those of Cyg X- 1 , 
iBelloni et alJ dl996l) found these to be very similar to the "in- 
termediate state" seen in these X-ray transients. In transients, 
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Fig. 10. Development of the power spectra near the flaring events seen in the ASM lightcurve. Their location in the ASM 
lightcurve is indicated in the insets. The panels show the data from a) 1998 July, b) 1999 December, c) 2000 October/November, 
and d) 2001 January/February. 



the "intermediate state" is typically seen during transitions be- 
tween the soft and the hard state, which fits the picture drawn 
above from our RXTE monitoring. 

Based on our observations of increased X-ray time lags 
during the transitional phases between the 1996 hard and soft 
state and during the ASM flares, we h ave previously called th e 
ASM flares "failed state transitions" JPottschm idt et al.l20 00). 
As outlined above, the behavior of the other timing quantities 
apart from the X-ray time lags now justifies our use of this 
term - during the "failed state transitions", the source changes 
its behavior compared to the normal hard state, by sometimes 
even entering the intermediate state, but a full soft state is gen- 
erally not reached, at least not on the sampling timescale of our 
RXTE monitoring. 

The availability of the different timing quantities for our 
monitoring observations also allows us to connect the increased 
X-ray time lag with the observed changes in the power spec- 
trum and the coherence function. As we mentioned above, 
the X-ray time lag var ies strongest in the 3.2-10Hz band 
( Pottschm idt et alj EoOO). The main contributions to the total 
rms variability in this frequency range are due to L2 and, to 
a lesser extent, due to L3. These two components are also the 
only broad noise components that are clearly present in the in- 
termediate states of 1999 December 05 and 2000 November 03 
power spectra (Fig.|9p and d), as well as during the tra nsition to 
the soft state in 199 6 (Fig. 2 o flPottschmidt et all200ll see also 
iBelloni et al.ll 996). Since the time lag behavior does not seem 



to reflect the transient behavior of the L3 component, and since 
the frequency range showing the enhanced lags is more consis- 
tent with that of L2, we tentatively identify the enhanced time 
lag during failed state transitions as being mainly due to L2. 
Since the soft state itself shows neither these broad noise com- 
ponents nor enhanced time lags, our identification provides a 
self-consistent picture. 

While most of the ASM flares are indeed "failed state tran- 
sitions", there is one flare during the time period covered here, 
in 2001 January, where the state transition did not fail, but 
where a soft state was reached for a brief time. The evolu- 
tion of the PSD shown in Fig.llOH displays again the vanishing 
of L\ and L4, with an "intermediate state" PSD observed on 
2001 January 29. From the point of view of our monitoring 
campaign, this behavior would have led us to classify this flare 
as yet another failed state transition. In addition to our RXTE 
monitoring, however, further pointed RXTE observations were 
performed in early 2001, that had a much bet ter sampling than 
our campaign. The analysis of these data bv lCui et alJ 1 2002) 
shows PSDs that are similar to those seen here, with the excep- 
tion of the data taken on 2001 January 28. Both the timing and 
spectral behavi or are very similar to the 1996 soft state, leading 
ICui etalJ l Eool to claim a possible short soft state of Cyg X- 1 . 

We note that in addition to the X-ray properties, the ra- 
dio flux was also peculiar in 2001 January. During the end of 
2001 January, the radio flux first dropped to about 30% of its 
typical value of ~15 mJy, before a radio flare was seen (Figs.|3^ 
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Fig. 11. Change of normalization of L\ as a function of en- 
ergy during the first half of 1998. Solid lines: Average of the 
normalization of L\, R\, for observations P30157/19 through 
28, before the 1998 July failed state transition. Dashed lines: 
Average of R\ for observations P30157/29 through 33, during 
the 1998 July event. 



and[4ji. In contrast to 1996, however, Cyg X-l was always de- 
tected in the radio, indicating that the system had not settled 
into the soft state as it did in 1996. In fact, our monitoring ob- 
servation of 2001 January 29 shows again the transitional PSD. 
As in 2001 November (Sect. 13.31 1. therefore, the switch from 
the soft state (if it can truly be labeled a soft state) back into the 
intermediate state occurred within one day. 

3.3. The 2002 Soft State 

Our RXTE observations performed since 2001 September give 
further credence to the identification of the presence of an "in- 
termediate state" during the transitions from the hard to the soft 
state. This time, from 2001 September until 2002 October, was 
the first time since 1996 that Cyg X-l exhibited again a full 
soft state behavior in the X-rays and in the radio emission. A 
detailed description of this s oft state will be th e subject of a 
later publication in this series JWilms etall2 003). here we will 
only briefly summarize the observations of the transition into 
the soft state during the end of 2001 as far as they are con- 
nected with the interpretation of the flares with the "intermedi- 
ate state" fSect.13.21. 

During 2001 September, the RXTE- ASM rate rose from its 
typical 20-30 cps to ~ 100 cps, with precursor flares in July and 
August (Fig.ll3>. Reaching its peak in mid-September, the soft 
X-ray flux first declined, and was found at a lower but still en- 
hanced ASM level of approximately 60 cps in 2001 December, 
the end of the time period considered here. However, in con- 
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Fig. 12. Evolution of the PSD during the monitoring observa- 
tions of 2000 November. The line style of the best fit compo- 
nents is the same as that in Fig. [5] 



trast to previous flares, this time a clear reduction of the radio 
emission was also observed, with almost zero flux for about 
one week in the beginning of October CFig.1131. After that time 
the radio emission slowly turned on again, although contin- 
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Fig. 13. RXTE ASM count rate and 15 GHz flux of Cyg X-l 
during the fall of 2001. The dashes denote the time of those 
pointed RXTE observations that were part of our monitoring 
campaign. 



ued flaring in both the radio and the X-rays was observed 5 . 
The combination of the strong ASM flux and the missing ra- 
dio emission is typical for a soft state. As shown in Fig. [H] 
pointed RXTE observations performed in 2001 October show 
a clear soft state behavior: comparing PSDs and time lags 
for the two monitoring observations of 2001 October 20 and 
2001 November 3 (dashes in Fig. 1 1 31 with typical data from 
failed state transitions and from the 1996 soft state shows that 
on 2001 October 20 the PSD had a clear / _1 shape, turn- 
ing over at ~10Hz, an X-ray time lag comparable with that 
measured during the standard hard state and an X-ray coher- 
ence that was close to unity. The time lag behavior, the unity 
coherence, and the PSD shape are typ ical for the soft state 
JCui et al.lll997tlPottschrnidt et alJl2000h . The X-ray spectrum 
of the 2001 October 20 PCA observation in the 3-20 keV range 
can be very well described by the empirical soft state spec- 
trum of ICui et"al] Jl997h . i.e., the sum of a black body and a 
broken power law. The broken power law has a photon index 
of Fi = 2.78 + 0.03 below, and T 2 = 2.10 + 0.02 above the 
break energy of E b = 1 1.1 ± 0.1 keV Of 2 = 27.5 for 31 de- 
grees of freedom, assuming an energy independent systematic 
uncertainty of 0.3%, a reasonable value for the current calibra- 
tion of the PCA). The temperature of the black body is mea- 
sured to kTsB = 0.36 ± 0.03 keV. The black body contributes 



5 We note that the radio monitoring has been more frequent than our 
pointed X-ray observations, and this is the first clear indication of the 
soft state in both the radio and the X-rays since the radio monitoring 
began. 
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Fig. 14. Time lags and PSD for the soft state of 2001 Oct com- 
pared with other typical soft state, hard state, and failed state 
transition data. 



20% to the total 2-10keV flux of 3.6 x 10~ 8 erg cirT 2 s. These 
spectral parame ters are similar to those of the 1996 soft state 
(Cui et al. 1997), although the observed 2-10keV flux is a fac- 
tor ~2 higher. We note that the 2001 October 20 data require an 
additional and very broad iron line at 6.4 keV (equivalent width 
~ 440 eV with a Gaussian width cr = 0.9 ± 0.1 keV) that had 
not been seen in 1996. 

On 2001 October 22, RXTE performed a second obser- 
vation close to this soft state observation. Here, the PSD 
slightly deviated from the soft state PSD, while the X-ray 
spectrum was still well described by the empirical soft state 
model. Shortly after this, during the next monitoring observa- 
tion on 2001 November 3, the PSD showed an "intermediate 
state" shape similar to that during the failed state transitions 
of 1999 November and 2000 December discussed above, i.e., 
with the L\ component missing and an increased power law 
component. Also pointing towards a transitional behavior are 
the increased time lag and the reduced X-ray coherence. Both 
the PSD and the X-ray time lag are similar to the behavior mea- 
sured on 2000 November 03, one year before the 200 1 soft state 
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(see Fig. 114-i . The X-ray spectrum is much harder and can be 
well described by the hard state spectral model. 

We conclude that in 2001 October and November, Cyg X-l 
was switching back and forth between the "classic" soft state 
and the "intermediate state". These changes from the soft state 
to the transitional behavior happened very quickly, within a 
few days at most. Furthermore, the very similar behavior of the 
source during the failed state transitions and the 2001 October 
transition into the soft state provides clear evidence that the 
ASM flares seen since 1996 are indeed "failed state transi- 
tions", as oppos ed to being brief " soft states". As we will show 
in a later paper dWilms et al.l2 003). after that episode Cyg X-l 
settled into a soft state until the end of 2002. 

3.4. Narrow Lorentzians 

In the previous sections we have described the overall evolu- 
tion of the four broad Lorentzians, which are responsible for 
most of the variability observed in Cyg X-l. Here we fin- 
ish our discussion by concentrating on the additional narrow 
Lorentzians, L a dd,i through L a dd,3, which were needed in about 
65% of all observations to completely describe the observed 
power spectra (Sect. 12.31 . As these components are always 
rather narrow, with Q » 2, we wil l follow the conventions out- 
lined, e.g., bv lBelloni et alJ d2002t) and call these components 
QPOs. As already mentioned in Sect. 12.31 these QPOs fall into 
two categories: rather narrow and clearly visible QPOs, L a dd,i 
and L a dd,2, with frequencies above roughly 1 Hz, and a broader 
structure, L a dd,3, at frequencies below 0.2 Hz. 

Our emphasis on modeling the broad Lorentzians led us to 
logarithmically rebin the measured PSDs, to the extent that a 
determination of the parameters of the narrow features was not 
always possible (see footnote[3}. For example, the logarithmic 
rebinning of the PSDs means that the QPO was concentrated 
in only a few frequency bins such that its Q-\alue could not be 
determined. Q-factors larger than 50 in Table[2]generaliy are an 
indication that these QPOs are concentrated in one frequency 
bin. The majority of the remaining 2-factors cluster between 
~2 and ~15, i.e., these QPOs are somewhat broader structures. 
Due to these limitations in our modeling, we can not quote 
formal uncertainties for the fit parameters of the QPOs. We 
will concentrate, therefore, in the following on the frequency 
behavior of the QPOs, keeping in mind the above problems 
with the detection and quantification of the QPO features. For 
consistency with the description of the broad Lorentzians, we 
will continue to use the peak frequency to characterize the fre- 
quency of the QPOs. Due to the large Q-values involved, the 
peak frequency is usually very close to the center frequency of 
the Lorentzian used to model the QPO. 

Turning first to the two narrow QPOs, Fig. ^] shows the 
distribution of the ratio between the peak frequencies of all nar- 
row Lorentzians to those of L\ and L^. Note that the majority of 
these additional PSD components are found between L\ and Li_. 
There are a few observations in which a QPO is present below 
the peak frequency of L\ , while no additional components were 
required above Li. As is clearly shown in Fig.ll5h. an apprecia- 
ble fraction of the QPOs are found at the first harmonic of the 



peak frequency of the L\ component, v\ (~ 20% of all detected 
QPOs have peak frequencies in the range 1.8 < v a dd/vi < 2.2). 
This enhanced fraction of QPOs at 2v\ suggests a relationship 
with the broad Lorentzian L\ . Taken at face value, one might 
interpret the QPOs in this frequency range as overtones of L\, 
and indeed L\ often shows a more complex structure which is 
modeled away by adding a Lorentzian to the model. We note 
that such a presence of possibly harmonically related c ompo- 
nents is not unexpected as it was seen before by Nowa kl (l2000h 
in his treatment of the PSDs of Cyg X-l and GX 339-4. 

We caution, however, that the "overtone" interpretation of 
this substructure is not the only interpretation. Since L\ is 
broad, one might expect its overtones to be broad as well, and 
this is not the case: the narrow QPOs with Q > 50, i.e., those 
QPOs which are only present in one frequency bin of our re- 
binned PSDs, have a distribution which is similar to that of 
QPOs with Q < 50 (Fig.ll5h. dashed lines). Furthermore, since 
for the peak frequencies of L\ and L2, v%lv\ ~ 8 (Table [3j, we 
have 2vi ~ V2/4, such that these QPOs could also be subhar- 
monics of Lq. We conclude, therefore, that the frequency dis- 
tribution of the QPOs is indicative of a relationship between at 
least part of these components and the broad Lorentzians, but 
that the quality of the current data is not good enough to distin- 
guish between the different interpretations. This is especially 
true given that the Lorentzian shape assumed for L\ through L4 
might only be an approximation to the real shape of the broad 
components that apparently comprise the PSD. 

Finally, we turn to the broad additional Lorentzian present 
at frequencies < 0.2 Hz, L a dd,3- In about half of the 17 observa- 
tions in which this component was present, an additional power 
law component was also required in the description of the data. 
Our data indicate a very loose correlation between the (J-value 
of i a dd,3 and its peak frequency. As shown in the previous sec- 
tions, PSDs containing a power law component are one indica- 
tor of (failed) state transitions. The presence of L a dd,3 in these 
PSDs thus indicates that the PSD during the failed state tran- 
sitions is slightly more complicated than the simple sum of a 
power law and the four broad Lorentzians and that the low fre- 
quency PSD is more stru ctured d uring these events (see, e.g., 
Miyamot o et alJl9 94 and Bellon fet alJl996l) . Finally, we note 
that we have not found statistically significant correlations be- 
tween the peak frequency of L a dd,3 and that of L\ through L4 or 
the narrow QPOs . The quality of our PSDs at low frequencies 
is not good enough, however, to allow a more detailed analysis. 

4. Discussion and Conclusion 

4.1. Summary 

The rapid high energy variability of galactic black holes is still 
not understood. The behavior of the PSD and the time lags dur- 
ing our campaign, however, gives some clue as to its origin. 

The most important result of this paper is that the shape of 
the PSD during the hard state and during the intermediate state 
can be productively interpreted with the decomposition of the 
PSD into the four broad peaked components, which we model 
as Lorentzians, and one power law. In contrast to earlier anal- 
yses, we were able to show that the decomposition holds for a 
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Fig. 15. a Number of detected thinner Lorentzian components as a function of v at jd/vi and b v a dd/v2- The dashed lines represent 
the same distribution for QPOs with Q > 50 only. The dotted line in panel a corresponds to a frequency ratio of 2; in panel b, 
dotted lines designate frequency ratios of 1:2 and 1:4 with respect to V2- 



very wide range of source parameters, even though the contri- 
butions of the power law and the Lorentzians are very different 
in the different states. Our result therefore adds further weight 
to the claim that the four broad Lorentzians are not just a con- 
venient description of the PSD ( merely more successful than 
the approximation of a broken power law), but that instead the 
broad components are the fundamental building blocks of the 
PSD of Cyg X-l and thus possibly also of the PSDs of other 
black hole candidates. We emphasize that a statistically signifi- 
cant discrimination between different shapes of the Lorentzians 
is not possibl e with the available data (see Sects. I2"31 and l3.4l 
and lBellonietall20 02). however, our study shows that any al- 
ternative building blocks for the PSD must at least be similar 
in shape to the broad components revealed by our analysis. 

Furthermore, there are two fundamental results from our 
campaign on the long term variability of the source. The first 
result is that the change of the general long term behavior of 
Cyg X-l from a "quiet hard state" to a "flaring hard state", with 
frequent intermittent failed state transitions after 1998 May, co- 
incided with a change in the PSD shape and amplitude. After 
1998 May, L3 was clearly much weaker than before. At the 
same time, the whole PSD shifted towards higher frequencies 
in a way that preserved the ratio between the peak frequen- 
cies, a larger fraction of PSDs showed low frequency noise, 
the X-ray lag increased, and the average X-ray spectrum soft- 
ened. The tendency for a softening of BHC spectra to corre- 
late with a frequency shift of charac teristic features o f the PSD 
has also been seen in other source s (di Matteo & Psaltis 1999; 
Gilfanov et all 1 19991 Eevnivtsev et alJ l200(A iKalemci et alJ 
2001: lNowak e t al. 2002), and thus might be considered an in- 
trinsic feature of all BHCs. 

The second result of our campaign concerns the change 
in the PSD shape during the X-ray flares. We interpret these 
flares as transitions from the hard state via the intermediate 
state into the soft state. Most of these transitions "failed", i.e., 
the transition stopped before the soft state was reached. That 
the "failed transitions" are caused by the same physical mech- 
anism as the normal transitions into the soft state was shown 



by the two instances during the campaign where a brief soft 
state was observed (Sect. 13. 2> . In these cases, before the soft 
state was reached and after it was left again, the source behav- 
ior was indistinguishable from the behavior during the "failed 
transitions". 

On the way into the intermediate state, L2 and L3 almost 
remained constant, while L\ - and, to a lower significance, also 
L4 - became remarkably weaker. At the same time, the X-ray 
spectrum softened and the whole PSD shifted to higher fre- 
quencies. Furthermore, the X-ray time lag increases in the 3.2- 
10 Hz frequency band during the transitions and the coherence 
function decreases. This is the frequency band where the PSD is 
dominated by the Li and, to a lesser extent, the L3 components. 
In those cases where the soft state was reached, the X-ray lag 
spectrum and the coherence function during the soft state were 
found to be very similar to that of the normal hard state, while 
the PSD, in contrast to the typical hard state, showed its char- 
acteristic f~ l shape with a cut-off at ~10Hz. During the soft 
state, no radio emission was detected. Such a detailed study of 
the hard state and transitions has to date only been possible for 
Cyg X-l. We note, however, that observations with a sparser 
cove rage show similar behavi or for GX 339-4 dNowak et all 
120021) and for XTE J 1650-50 dKalemci et al.l2003l) . 



4.2. Relationship between Neutron Star and Black 
Hole Sources 



Apart from the analysis of Nowak (2000), Lorentzian mod- 
eling has also been used for describing the power spectra of 
neutron star sources. A recent summary of Lorentzian fits to 
data from neutro n stars and bla ck holes has been recently 
published by Bell oni et alJ ( 120021) . These authors base their 
work on their analysis of RXTE observations of six neutron 
star sources (IE 1724-3045, SLX 1735-269, GS 1826-24, 
XTE Jl 118+480, Cir X-l, and GX 17+2), as well as on 
published results fr om the black holes Cyg X-l, GX 339- 4 
dNowak et al.l2002l) . and XTE J 1550-05 dHoman et al.l200ll) 
and the neutron star sources 4U 1915-05 (Boirinetal. 
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Fig. 16. Correlations b etwee n the peak frequencies of the broad noise components, L\-L\ and the frequency correlations dis- 
cussed by Psal tis et 101 (17999. solid lines are v v = 42Hz(v r /500Hz) 095±016 ). a) v\ versus V2, b) V2 versus V3, c) v\ versus vj. d) 
v\ (crosses), V2 (diamonds), and V3 (circles) versus the highest peak frequency, V4. Only data taken in P40099 and P50110 are 
shown in panel d to avoid biasing due to the 32 Hz frequency cutoff in the P30157 data. 



200 0| and 4U 1728-34 JPi Salvo et alJ 1200 ll) . Furthermore, 
van Straaten et alJ J2002I) have recently published an extensive 
an alysis of data from 4U 06 14+09 and 4U 1728-34. 

iBelloni et ail J2002I) and lvan Straaten et alJ J2002I) find that 
one to six Lorentzians are required to obtain a good de- 
scription of the PSDs of these sources. A slight complication 
when com paring these results t o ours is that lvan Straaten et alJ 
( 2002) and Be lloni et alJ ( 120021) mainly use zero frequency cen- 
tered Lorentzians, i.e., shot-noise profiles, and not Lorentzians 
with non-zero center frequencies. As we describe in Sect. 12.31 
our choice appears to result in a slightly lower y 2 (but note 
the discussion in Sect. 2 of IBelloni et alJl2002t) . Despite the 
different fit functions, however, the behavior of the maxi- 
mum frequenci es in both cases is still c omparable. In terms 
of ascending v, Ivan Straaten et alJ J2002I) call L\ and L2 the 
band-limited noise and zero-centered Lorentzian, L3 the low- 
frequency Lorentzian, an d L4 th e hectohertz Lorentzian. In the 
notation of Bellon iet alJ ll2002l) . L\ corresponds to Lb as this 
is the Lorentzian describing the break frequency of the PSD. 
The Lorentzians L3 and L4, which describe the lower and up- 
per characteristic frequencies in the tail of the PSD, are called 
L\ and L u . In many neutr on star PSDs t he Lq component has a 
rather large Q, therefore, Belloni et al. (2002) identify L2 with 
the low frequency QPO, Llf- 

Note that the for all sources the "correct" identification 
of these components is difficult. For example, L\ and Li are 



not present together in all data sets such that one might mis- 
int erpret the band limited noise (to stay in the terminology 
of Ivan Straaten et alJ Eool as the zero-centered Lorentzian. 
Similar problems also apply for Cyg X-l: had the spacing 
of our observations during the (failed) state transitions been 
coarser, we would not have been able to correctly identify the 
Lorentzians during the state transitions. 

At higher frequencies, additional narrow Lorentzians are 
required to describe the kilohertz Lorentzians in neutron star 
sourc es, however, the ir connection to the PSDs of BHCs is less 
clear. Nowak (2000) argue tentatively that L4 could be iden- 
tified with the upper kilohertz QPO , however, the obs erved 
behavior of this component bv lvan Straaten etafl J2002I) casts 
some doubt on this interpretation: while neutron star kHz QPOs 
are varying in frequency, L4 clearly doe s not (e.g., Fig. (6j. On 
the other hand. Ivan Straaten e t al. (2002) point out that there is 
a broad feature around 100 Hz, the "100 Hz bump", which is 
constant in frequency and therefore a more probable counter- 
part to L4. 

There are two main correlations among frequency com- 
ponents that have been claimed for neutron star and black 

I 1 1 I 

hole sources (see Belloni et al. 2002 for a recent summary). 



Wiinands & van der Klis ( 1999) have pointed out a correlation 
between the PSD break frequency and the low frequency QPO 
of X-ray binaries. These two frequencies might be related to vi 
and V2 in our classification scheme. Psalti s et al.l 11999) have 
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pointed out a correlation between the low frequency QPO and 
the lower frequency kHz QPO, i.e., frequencies that could be 
identified with V2 and V3 in our scheme. The latter correla- 
tion especially relies on identifying these features in multiple 
sources. The origin of these correlations is unknown. 

In Fig. we plot various combinations of one frequency 
versus another, and show that with the exception of v\ all 
frequency combinations show some sort of correlation (this 
is just a consequence of the constancy of the frequency ra- 
tios shown in Fig. Similar correlations for several sources 
are also shown by Bell oni et alJ 12002. Figs. 11 and 12), who 
point out that - while all different sources seem to show cor- 
relations between individual Lorentzians - these do not seem 
to completely agree from source to source. Given the above 
two correlations, t he one that f its our data best is the correla- 
tion described by IPsaltis et al.l dl999h . albeit applied to the v\ 
versus V2 correlation. We note, however, that in this particu- 
lar range of frequencies, the vi versus vi correlation cited by 
Wiinands & van der Klis (199 9J ) is nu merically not very differ- 
ent than the V2-V3 correlation of Psalti s et al.l I I1999I) . Likewise, 
the V2-V3 data presented here are marg inally consistent with the 
correlation discussed bv lPsaltis et alJ ( H"999l) . Although the in- 
dividual components discussed here span more than a decade 
in frequency, this is not large enough a range to identify unam- 
biguously given features with the t hree decade wide c orrela- 
tions hypothesized and discussed by Psalti s et alJ II19991) . Still, 
correlations among the frequency components clearly do exist, 
and point to possible physical scenarios, as we further discuss 
below. 

4.3. A Possible Physical Scenario 

We now discuss the empirical evidence from our campaign 
in terms of more physical models for the accretion pro- 
cess in a black hole system. The hard state X-ray spectrum 
is thought to be caused, at least in part, by soft photons 

from the accretion disk which are Compton upscattered in a 
1 11 1 1 1 1 || 1 

hot electron gas (Sunvaev & Trumpei 1979; Dove et al. 1998; 

IPoutanen & Svenssonlll999l and references therein). In this 
picture, state transitions are caused by the disappearance of 
this Comptonizing medium. The cause for this disappearance 
is unknown, however, due to the luminosity difference be- 
tween the hard and the soft state, many workers have as- 
sumed that the Compton corona is produced by some physi- 
cal process which can work only at mass accretion rates be- 
low a ce rtain t hreshold rate such as the magnetorotation al in - 
stability ( Balb us & Hawlevll99ll) . ad yection fesin et al.1 19971) 
or some kind of coronal ou tflow {Merloni & Fabianl |2002; 
iBlandford & Begelmanl ll999. and references therein), while 
others have claimed the presence of an up to now un- 
known additional parameter which triggers the state transi- 
tions faoman et all200ll) . A similar picture of changing coro- 
nal properties should also hold true within the hard state it- 
self - although the Compton corona exists throughout the hard 
state, phases where a different X-ray spectrum and PSD are 
observed mig ht well correspond to differen t coronal geome- 
tries (Nayakshin et al. 2000; iLiu et aljfl999l lEsin et al.Ml997t 



ISmifh et aljEooi e.g. ,). Such a change in the coronal config- 
uration could explain the reduction of the L3 component after 
1998 May, as is also evidenced by an overall softening of the 
spectrum since that time. 

We note that forCyg X-l the luminosity differe nce between 
the ha rd and soft state is rather small (~35%; IZhang et alJ 
Il997ah . Furthermore, a power law component is always seen 
during the soft state. This suggests that the soft state in Cyg X- 
1 might not be a "normal" soft state as that seen in other black 
holes, but rather closer to the very high state than the canon- 
ical soft state is. Furthermore, there are cases where the soft 
state has not a higher luminosity but a lower luminosity than 
the hard state. For GRS 1758-258, an intermediate state has 
even been observed at a higher luminosity than the hard state, 
whil e the soft state a ppears to be at a much lower luminos- 
ity JSmifh et all200ll) . A possible explanation for these devia- 
tions from the canonical picture of the soft state having a much 
high er luminosity than the h ard state might be an hysteresis ef- 
fect (Miyam oto et aljfl99 5). however, we note that regardless 
of the physical cause for deviations the observational evidence 
still points towards changes in the accretion disk geometry. 

That the accretion disk geometry changed around 
1998 May might also be the cause for the shift of the peak 
frequencies towards higher frequencies observed during that 
time. Even in prior analyses where PSDs of BHC have been 
described in terms of power law fits and "shot noise models" 
jLochnereta"il l991. and references therein), it has been pre- 
sumed that the characteristic fre quencies can be attribu ted to 
time scales of the accretion disk. IChurazov et all J200ll) have 
argued that PSD "breaks" are related to the size scale at which 
the accretion flow transits from a thin disk into a geomet- 
rically thick, hot inner corona. Di Matteo & Psaltis (1999) 
have argued that the dynamic range over which such a tran- 
sition varies is somewhat limited since the observed frequen- 
cies themselves, although variable and correlated with spec- 
tral hardness, only span a limited rang e. Since flow fre quencies 
scale as 5v/v ~ (H/R) 2 (e.g. Jdi Matteo & Psaltis! 19991 and ref- 
erences therein), if the Lorentzians are associated with resonant 
effects within the accretion disk, they must originate from a re- 
gion that is geometrically thicker than the thin accretion disk, 
i.e., a comparably extended region such as the accretion disk 
corona. 

A simplified analysis of the characteristic frequen- 
cies of such a transition region has been presented by 
Psal tis & Normanl J200 lh . These authors note that if the outer 
disk acts as a "noise source", the transition region acts as a 
"filter", and the inner corona acts as a "response" which pro- 
cesses the filtered noise into observable X-ray variability, then 
one expects to observe a number of quasi-periodic features. 
These features would be related to radial, vertical, acoustic, 
rotational, etc., oscillations of the transition region. Many of 

—3/2 

these frequencies would be expected to scale as R T , w here 
Rt is t he transition region radius. More recently. iNowak et alJ 
(2002) have shown that such a scaling might exist for hard 
state observations of GX 339-4. Specifically, they show that 
in "sphere and disk" coronal models, the "coronal compact- 
ness" - where higher compactness, f c , means harder spectra 
- scales as l z rc R T . Furthermore, the characteristic PSD fre- 
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quencies, measured in a manner similar to that used to describe 
the Cyg X-l PSDs presented here, scale as £~ 3/2 . Both Cyg X-l 
and GX 339-4, in their hard states, show that spectral softening 
is correlated with higher PSD frequencies, although we have 
yet to characterize th e Cyg X- 1 spectra in terms of the coronal 
models presented by Nowak et al. (2002). Such an analysis is 
currently underway (Gleissner et al. 2003). 

Further hints towards attributing at least some of the 
Lorentzians to the Compton corona come from the PSD be- 
havior during the transitions into the intermediate or the soft 
state. We hypothesize that the variability causing L\ and possi- 
ble L\ is also associated with the accretion disk corona, since 
these components - like the corona - appear to vanish during 
the transition. Since hi and L3 remain present in the 3.2-10 Hz 
band, we are tempted to attribute the increased X-ray time lag 
in this f requency band to these components. As discussed by 
iNowakl §000), the overall measured time lag may in reality 
be a composite of time lags inherent to each individual PSD 
component. Whereas the composite time lag might be fairly 
small, the time lags associated with the individual PSD com- 
ponents might be somewhat larger (specifically, see Fig. 5 of 
iNowaktZO OO). The fairly long intrinsic time lags of these indi- 
vidual PSD components would only be revealed in the absence 
of the other PSD components. Such a picture might explain the 
larger X-ray time lag during state transitions: Here, the Li and 
L3 components, which sit in the 3.2-10 Hz band, contribute a 
larger relative fraction of the total source variability. If their 
intrinsic lag is larger, a larger lag would be observable during 
state transitions when the contribution of L\ and La is reduced. 
We note that such a picture of a "composite X-ray time lag" is 
somewhat suggested by the observational shape of the time lag 
spectrum, where "jumps" between different frequency bands of 
roughly constant time lag seem to correspond to the overlap re- 
gions of the broad Lorentzian s JNowak et alJll999bt IPoutaner] 
l2002HPottschmidt et all200(il) . 

While the origin of the time lag is not yet understood, it is 
possible that the time lag is somehow related to a character- 
istic length scale of the medium which produces the observed 
photons. If this is true, then the larger lags during state transi- 
tions may be showing that some part of the source geometry 
changes and becomes larger during the transitions. It is im- 
portant to note that although there is no radio emission dur- 
ing the soft state, there does seem to be a slight increase of 

the radio flux at least during some of successful transitions 

1 11 1 1 I 

into and out of the soft state (Corbel et al. 2000; Zhang et al. 

Il997bl) . Furthermore, intermittent, rapid radio flaring activity 

also seems to be associated with the (failed) transitions (Pooley, 

priv. comm.). 

In general, we see that the time lags are correlated with 
the PSD frequencies and correlated with the photon index, T 
(Fig.|7Ji in Sect.EJ. F or GX 339-4, spectral fits with the coro- 
nal model of Coppi ( 1999) - where a spherical corona with 
seed photons distributed throughout it according to the diffu- 
sion equation is modeled - suggest an inc reasing radius with 
a softening spectrum jNowak et alJ 120021) . cons istent with a 
geom etric interpretation of increased time lags. iNowak et alJ 
(2002) therefore argued for a jet-like model for the hard state 
of GX 339-4 where the radius of the base of the jet followed 



the radius-frequency correlations suggested by the "sphere and 
disk" coronal model (i.e., smaller base radius yielding softer 
spectra), while the vertical extent of the corona/jet followed the 
correlations suggested by the Coppi coronal model (i.e., more 
vertically extended coronae yielding softer spectra). Thus, fre- 
quencies might be generated at the base of the jet, while time 
lags could be generated via propagation along the jet. Given 
the similar, but much more detailed, correlations observed here 
in Cyg X-l it is tempting to ascribe the same type of model 
to these data. Whether the analogy to GX 339-4 holds in de- 
tail w ill depend upon the results of detailed coronal model fits 
JGleissner et all2003h . 

Along these lines, b ase d mainly on observat ions of micro- 
quasars, Fended (120011) and tFender et al.1 (11999) have recently 
attempted to unify the radio and X-ray observations into a ge- 
ometrical picture of the region surrounding the central black 
hole. We have previously suggested that the enhanced X-ray 
lags during the "flare s" observed in Cyg X-l s eem to add cred- 
ibility to this picture dPottschmidt et all200ll) . During the nor- 
mal hard state, there is a weak outflow of material from the 
region surrounding the black hole. This outflow is responsi- 
ble for producing the observed synchrotron radio emission. We 
think it likely that the base of this outflow coincides with the ac- 
cretion disk "corona" that produces the hard state X-ray spec- 
trum and is, given the results above, also responsible for part 
of the observed X-ray variability. During the "flares", part of 
the corona gets ejected. Thus the mass of the radio outflow is 
temporarily enhanced, increasing the radio luminosity. As the 
corona is disturbed, whatever resonant mechanism producing 
the L\ and possibly also the La variability will be disturbed, 
resulting in a decrease of its contribution to the total rms vari- 
ability. Furthermore, if this "corona" is "stretched" because of 
the ejection, it seems likely that a temporary increase in the 
X-ray lags should be associated with the outflow as the charac- 
teristic length scales increase. During the soft state, the corona 
has fully vanished such that no radio outflow is observed and 
the X-ray spectrum is soft. If L2 and L3 are produced at the 
corona-disk interface, they also vanish. 

This very rough picture leaves many detailed questions 
open, such as the explanation for the strong domination of the 
/ _1 power law during the soft state, its cutoff frequency, the 
generation of the intrinsic lag of the Lorentzians, or the phys- 
ical mechanisms responsible for the PSD itself. However, we 
believe that monitoring campaigns such as this one, that are 
able to assess the long term, systematic changes in accreting 
systems, will continue to constrain the parameters of the ob- 
served physical systems and will enable us to switch from de- 
scribing the "accretion disk weather" towards understanding 
the "accretion disk climate". Equally importantly, X-ray as- 
tronomy is entering a phase of high spectral resolution obser- 
vations with such instruments as Chandra and XMM-Newton. 
The nature of these instruments does not allow detailed and in- 
tensive monitoring campaigns such as this one. However, given 
the continued operation of the RXTE-ASM and judicious use 
of simultaneous pointed RXTE observations, campaigns such 
as this will allow these individual high resolution observations 
to be placed within their proper 'global' context. 
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